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SEMI-ANNUAL  TECHNICAL  SUMMARY  REPORT 
June  1966  -  December  1966 

1.  INTRODUCTION  ■ 


A.  Extension  of  Contract. 

Under  a  Supplemental  Agreement,  effective  1  January  19^7, 
Contract  AF  49  (638)-1506  has  been  extended  until  31  December 
1968. 


Hence,  this  report  is  not  a  Final  Report  of  research  but 
the  Semi-Annual  Technical  Summary  Report.  Nevertheless,  we 
thought  it  appropriate  to  bring  together  here  in  summary 
form  the  investigations  undertaken  during  the  first  two  years 
(1965  -1966)  of  the  contract. 


B.  Research  Scope.  ^  ^  .  j  ■  |  f  ■  i  >  <  r  ^ , 

The  contract  is  entitled  ’’Research  on  the  Central  California 
Large-Scale  Seismic  Array".  The  objectives  of  the  study  are  as 
follows : 

1.  Maintain  and  operate  the  telemetered  network  of  stations 
in  the  seismically  active  Coast  Range. 

Jp/.  4 Conduct  research  on  the  seismic  data  recorded  from  this 
network.  This  research  should  include  the  following: 


a.  Focal  mechanism  and  aftershock  characteristics  of 
small  earthquakes 

b.  Propagation  of  S  waves  across  the  network 

c.  P-wave  coherence  from  teleseisms 

d.  Improved  travel -time  curves.  ;  } 


2.  SUMMARY  OP  RESEARCH  ACCOMPLISHED 


A  telemetry  network  of  ten  telemeter  stations  has  been  kept  in 
operation  in  1965-1966-  The  arraj  is  240  km  long  (Priest  Valley  to 
Berkeley)  and  l60  km  wide  (Granite  Creek  to  Jamestown).  This  is  the 
heart  of  16  seismographic  stations  operated  by  .he  University  of 
California,  Berkeley. 

The  network  has  allowed  high-quality  simultaneous  recording  of 
seismic  signals  on  film  and  magnetic  tape.  An  extensive  library  of 
recorded  wave  trains  from  earthquakes  (local  and  distant)  and  explosions 
has  been  built  up. 

A  careful  study  of  the  morphology  of  earthquakes  in  central  and 
northern  California  has  been  carried  out;  the  telemetry  facility  allows 
locations  tc  be  determined  with  more  precision  than  previously.  Atten¬ 
tion  has  been  given  to  the  variation  of  low-energy  (magnitudes  between 
2.5  and  5)  earthquakes  with  time  and  in  relation  to  mapped  fault  traces. 

Three  earthquake  sequences  associated  with  main  shocks  of  magnitude 
5  and  greater  have  been  analyzed.  These  and  earlier  sequences  begin  to 
suggest  a  division  into  two  types  of  behavioral  patterns. 

Considerable  work  has  been  done  on  the  propagation  of  P,  pP  and 
FKP  waves  from  teleseisms  across  the  array.  We  were  among  the  first  to 
report  the  existence  of  large-scale  azimuth  au.  mly  functions  for  the 
travel -times  of  P  waves  between  source  and  station.  Anomaly  functions 
for  most  of  the  Berkeley  stations  have  been  measured;  the  geological 
consequences  are  being  workt i  out. 

Again,  we  were  among  the  first  to  estimate  (and  publish)  mean  azi¬ 
muth  (and  slowness)  anomaly  functions  for  a  seismic  array.  The  functions 
have  been  incorporated  into  the  routine  rapid  location  of  teleseisms  car¬ 
ried  out  at  Berkeley.  We  have  found  that,  for  example,  surprisingly  close 
locations  for  Longshot,  the  French  nuclear  tests  near  Tahiti  and  the 
U. S. S.R.  nuclear  tests,  can  be  computed  with  the  Berkeley  network  first- 
arrivals  if  the  anomaly  functions  are  taken  into  account. 


. .  ■iHtUrtwaig  i M-  LB— J—  j^ilhklto^WW»*«WMiai  -  ■  -  -mmirnm  (teillltiii!  aHHW 
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A  comparison  between  the  signatures  of  the  Chase  V  explosion  and 
an  earthquake  near  Chico,  California  has  drawn  attention  to  a  number  of 
different  elastic  wave  properties.  These  events  were  of  similar  magni¬ 
tude  and  were  recorded  on  the  augmented  network  of  seismogiaphs  specially- 
planned  by  us  for  the  Chase  V  blast. 

Our  research  is  an  directly  or  indirectly  germane  to  the  problem 
of  the  nature  of  the  crust  and  upper  mantle  in  this  part  of  the  globe. 

We  have  used  travel-times  of  P  and  S  from  quarry  blasts,  local  earth¬ 
quakes,  Nevada  atomic  explosions  and  Chase  V  to  determine  the  local 
crustal  structure.  We  have  extended  the  traditional  analysis  to  include 
the  use  of  S  particle  motion,  crustal  transfer  function  for  P  waves  of 
high  frequency,  and  the  station  azimuth  anomaly  functions. 
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4.  NETWORK  OPERATION 
W.C.  Marlon  and  R.W.  Sell 

The  Berkeley  network  of  stations  (telemetered  and  photographic)  has 
continued  to  operate,  with  preventive  and  corrective  maintenance,  in  a  sat¬ 
isfactory  manner..  Down  time  for  the  last  two  years  has  been,  on  the  average, 
less  than  one  per  cent  per  station;  all  records  ore  in  archive  storage. 

Table  1  lists  the  Berkeley  net  as  of  December  31>  1966* 

Magnet ic  Tape  Rec ording .  The  magnetic  tape  data  acquisition  and  anal¬ 
ysis  system  was  brought  up  to  full  capacity  in  June  of  1965  by  the  addition 
of  nine  recording  and  playback  channels  of  electronics,  a  time  code  genera¬ 
tor,  generating  B.C.D.  code  in  I.R.I.G.C.  format,  a  time  code  translator 
and  tape  search  unit,  and  filters  for  data  processing  with  visual  playout. 

The  translator  will  decode  Vela  Uniform  as  well  as  I.R.  I.G.C.  codes.  Table 
2  shows  station  tape  channels  by  channel  number  and  date,  listing  all 

cnanges  during  the  two-year  period. 

Seismic  events  recorded  on  tapes  are  dubbed  onto  the  array  library 
tapes,  which  now  contain  a  large  collection  of  local  earthquakes  and 
teleseisms. 

A  device  was  constructed  (design  assistance  was  given  by  Dr.  P.W. 

Rodgers)  to  resolve  the  outputs  of  the  horizontal  components  of  the  analog 
magnetic  seismograms  into  longitudinal  and  transverse  components  (Figure  l). 
The  system  utilizes  three  operational  amplifiers  as  analog  multipliers  and 
integrators  to  present  true  ground  particle  displacement  for  direct  conver¬ 
sion  to  motion  diagrams.  J.R.  Filson  has  used  this  device  successfully  for 
detailed  studies  of  the  particle  motion  of  the  S  phase  from  local  earthquakes. 

Changes  in  Central  California  Seismic  Array.  As  of  January  1,  1965  the 
following  telemeter  array  stations  were  in  operation: 


VIN  -  Vineyard 

SCC  -  Santa  Cruz 

PRI  -  Priest 

PR3  -  Paraiso 

CNC  -  Concord 

LLA  -  Llanada 

MHC  -  Mount  Hamilton 

$ 
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On  February  4,  1965  a  visible  recording  seismograph  was  installed  at 
the  U.S.  Coast  and  Geodetic  Survey  station  at  Ukiah: 

UKIs  39°  08.2'  N,  123"  12.6'  W,  Elevation  199  nu 

A  14  kg  Benioff  vertical  seismometer  amplified  by  a  photo-tube  ampli¬ 
fier  with  °,2  sec  galvanometer  and  filtered  by  a  special  low-pass  RC  filter 
designed  by  P.W.  Rodgers  was  recorded  on  a  helicorder  (Figure  2). 

A  new  telemetered  station  known  as  Granite  Creek  was  installed  February 

2 6,  1965: 

GCC:  37°  01.8'  N,  121°  59.8*  W,  Elevation  122  m. 

This  location  is  about  2  km  north  of  the  Santa  Cruz  (SCC)  station.  After 
this  location  had  greatly  improved  ~ignal/noise  characteristics,  the  Santa 
Cruz  station  was  discontinued  March  3>  1965- 

A  three -component  set  of  matched  Sprengnether  seismographs  (0.46  kg, 

Tq  =  2  sec)  were  installed  after  careful  calibration  at  the  Vineyard  tele¬ 
meter  station  March  24,  1965.  Signals  were  recorded  on  magnetic  tape  to 
provide  data  for  the  study  of  S  waves  in  the  crust  (Figure  3). 

A  new  telemetered  station  was  installed  on  the  San  Francisco  Peninsula, 
about  20  miles  south  of  that  city,  on  April  29,  19&5.  The  new  station, 
known  as  Pilarcitos  Creek,  will  provide  the  necessary  control  for  epicenter 
determination  in  the  Bay  Area: 

PCC:  37°  30. O'  N,  122°  22.9'  W,  Elevation  91  m. 

On  July  1,  1965  the  two  horizontal  Sprengnether  seismometers  at 
Vineyard  were  replaced  with  Willmore  seismometers  (4.75  Eg>  Tq  =  3.0  sec) 
(Figures  4  and  5). 

A  new  telemetered  station  was  installed  at  the  Jamestown  seismographic 
station  September  16,  1965 : 

JAS:  37°  56.8'  N,  120°  26.3’  W,  Elevation  457  m. 

The  three -component  net  of  seismometers  at  Vineyerd  w».s  removed  on 
October  25,  1965  and  the  station  was  put  on  standard  'operation  with  a  verti¬ 
cal  14  kg  Benioff  seismometer.  Due  to  construction  at  the  winery  the  Vineyard 
station  was  moved  to  a  temporary  site  on  the  Harris  Ranch  on  March  17,  19&5: 

HRC:  .36°  46.2'  N,  121®  24.8'  W,  Elevation  228  m. 
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Upon  completion  of  the  underground  vault  at  the  San.  Andreas  Geophysical 
Observatory  the  station  was  moved  again  on  July  11,  1966: 

SAO:  36°  ^5.9*  N,  121"  26.7*  W,  Elevation  230  m. 

A  raicroearthquake  channel  with  a  6-20  cps  bandpass  and  a  helicorder 
magnification  of  8  million  at  20  cps  started  telemetering  from  SAO  on 
September  l6,  1966.  Near  events  below  the  detection  threshold  of  the 
standard  system  are  being  clearly  recorded  on  the  microearthquake  channel 
at  a  rate  of  several  per  day  (Figure  6). 

The  visible  recording  seismograph  was  removed  from  Ukiah  November  6, 
1966.  The  Concord  (CNC)  signal  telemetered  to  Berkeley  was  discontinued 
effective  December  7>  1966.  Recording  is  now  being  done  locally  on  a 
helicorder. 

Mobile  Operations.  An  earthquake  of  magnitude  3.7  occurred  on  June  28, 
1965  in  the  vicinity  of  Calaveras  Reservoir,  latitude  37°  31 1  N,  longitude 
121°  45'  W.  Mobile  equipment  was  sent  out  and  established  on  location  a 
few  hours  after  the  event. 

An  'orthquake  of  Richter  magnitude  clobe  to  5  occurred  on  September 
10,  1965  in  the  vicinity  of  Mt.  Diablo,  about  25  km  east  of  Berkeley.  One 
foreshock  and  a  sequence  of  about  80  aftershocks  were  observed.  A  mobile 
station  was  operated  for  10  days  in  the  epicentral  area. 

An  earthquake  of  Richter  magnitude  3.6  occurred  on  July  18,  1965  in 
the  area  southwest  of  Vineyard.  The  mobile  unit  was  sent  out  but  no  sig¬ 
nificant  aftershock  sequence  was  recorded. 

An  earthquake  of  Richter  magnitude  4.1  occurred  on  January  17,  1966 
near  San  Felipe  Lake.  A  recording  truck  was  sent  out  January  l4  and 
recorded  the  main  shock  and  several  dozen  aftershocks. 

The  mobile  unit  was  operated  at  the  former  U.C.  seismographic  site 
at  Shasta  for  the  Chase  V  explosion  on  May  24,  1966, 

The  mobile  station  operated  in  the  epicentral  area  of  the  Parkfield 
earthquake  of  June  28,  1966  tnrough  July  9>  1966.  It  was  also  used  to 
record  recent  calibration  blasts  by  the  U. S.G.S.  near  Parkfield. 

The  mobile  station  was  sent  to  the  area  of  Ires  Pinos,  7  miles  south¬ 
east  of  Hollister  to  record  the  NTS  blast  of  December  20,  1966  at  the 
request  of  the  U. S.G.S.  of  Menlo  Park. 
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Special  Projects.  An  additional  telemetry  channel  with  0-20  cps 
response  was  operated  from  Pilarcitos  (ICC)  from  December  30,  1965  to 
January  29,  1966-  The  output  from  a  Wi.’lmore  1  cps  seismometer  was  ampli¬ 
fied  by  an  operational  amplifier  and  transmitted  on  this  channel.  Late  in 
January  a  swarm  of  micioearthquakes  (l.  5  sec  S-P  time)  were  recorded. 

On  February  4  a  1  km  triangular  array  of  similar  instruments  was  instaxied. 
By  March  12,  when  the  array  was  re:  oved,  no  more  such  earthquakes  had  been 
observed. 

An  attempt  to  record  the  very  low  velocity  (l  km /sec)  surface  waves 
across  the  Great  Valley  to  Jamestown  (JAS)  from  quarry  blasts  on  the  west 
side  was  made.  A  4-sec  Willmore  seismometer  was  installed  with  a  DC 
transistorized  amplifier  at  JAS  from  Fevruary  22,  1966  to  March  15,  1966- 
Since  the  surface  waves  were  in  the  5*8  sec  microseism  range,  it  was  found 
that  the  largest  quarry  blasts  generated  insufficient  amplitude  in  the  sur¬ 
face  wave  train  for  good  detection  in  the  Sierra  foothills.  In  fact,  the 
study  showed  at  least  a  magnitude  3  earthquake  is  required  as  a  source  of 
perceptible  waves  at  JAS. 

An  additional  telemetry  channel  was  operated  from  JAS  on  December  20, 
1966.  The  output  from  a  Press-Eving  vertical  seismometer  set  at  10  sec 
free  period  was  amplified  by  an  operational  amplifier  and  transmitted  on 
this  channel.  The  system  was  operated  at  very  low  magnification  for  the 
unlerground  explosion  on  that  date,  at  the  Nevada  Testing  Site. 
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Table  1 


SEISMOGRAFHIC  STATIONS  OPERATED  BY  OR  AFFILIATED  WITH  THE 
UNIVERSITY  OF  CALIFORNIA 

December  3lj  1966  . 


Instrument  Types 

B  =  Benioff  100  KG 

b  =  Benioff  14  KG 

P  =  UED  Long-Pe.  iod 

S  a  Sprengnether  Long-Period 

G  =  Geotech  Long-Period 

W  =  Wcod-Anderson 

L  =  Loucks-Omcri 

s  =  Sprengnether  Short -Period 

Z,N,E  =  Components  operating.  Vertical, 
North -South,  East -West 
nv,  ne  =  Components  operating, 

NU5°W,  N450E 

H  =  6-20  cps  Mic r oearthquake  system 


Recording  Method 

S  ■  Smoked  paper 

P  ■  Photographic 

H  r>  Helic  order 

Hm  ■  Helicorder  telemet  y  monitor 

sequenced  daily 
V  ■  Other  visible  'Mary; 

F  ■  l6  mm  film  (D-'  ilocorde.r) 

T  «  Magnetic  tape 

Magnllication 

(V)  ■  l6  mm  film  vievjd  at  20X 


Location 

Instruments 

Magnification 

Sta¬ 

Area 

Elev. 

Latitude 

Longitude 

Recording 

tion 

(Calif. ) 

(m) 

(N) 

(w) 

Type 

Ts 

Ts 

Max.  1  cps 

Method 

ARC 

Areata 

59 

40'  52.’6 

124°  04.' 5 

bZ 

1.0 

0.2 

5900 

P 

(.877) 

(.075) 

WNE 

0.8 

- 

2800 

P 

BRK 

Berkeley 

81 

37°  57.’ 4 

122°  15:6 

WKE 

0.8 

- 

^00 

P 

(E.  S.  B. ) 

(.873) 

(,260) 

WNE 

0.8 

- 

4 

P 

BZ 

1.0 

0.2 

24000(V) 

F,  Hm 

2400(V) 

(Hav.  H. ) 

BZ 

1.0 

8.0 

3000+ 

V 

.,,.nv 

PZ„ 

ne 

30 

- 

200  & 

20  v/cm/ 

T,  H(Z) 

sec 

.04-10  cps 

PZ 

15 

30 

1000 

H 

Pne 

45 

300 

650 

(40  sec) 

P 

BKS 

Berkeley 

276 

37°  52.’6 

1226  l4.'l 

WNE 

0.8 

2800 

P 

(Strawberry 

(.877) 

(.235) 

BZNE 

1.0  0.75 

35000  25000 

P 

Canyon) 

SZNE 

15 

100 

3000 
(15  sec) 

P 

CJC 

Concord 

36 

37°  58:1 

122°  04.' 3 

BZ. 

1.0 

0.2 

40000 

H 

(.968) 

(.072) 

Tabl«  1  (continued) 
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Location _  Instruments  Magnification 


Sta¬ 

tion 

Area 
(Calif. ) 

Elev. 

(m) 

Latitude 

-(*) 

Long itude 
(W) 

Type  Ts 

T  Max. 

g 

1  cps  Recordin« 
p  Method 

FRE 

Fresno 

88 

36°  46.'0 

119°  47.'8 

sZHE 

2.0 

2.0 

6500(ra^ 

p 

(.767) 

(.797) 

109G0(Z) 

JCC 
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122 

37°  01.'8 

121°  59.'8 

bZ 

1.0 
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4oooo(v) 

F,  Hm 

Creek 

(.030) 

(.996) 
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1.0 

0.75 

260000 
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(.438) 

bZ 

1.0 

0.2 

6coooo(v) 

T,F.K 
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Llanada 

475 

36°  37  .'0 

120°  65:6 

bZ 

1.0 

0.2 
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F,  Hm 

(.617) 

(.943) 

MHC 

Mt. 

1282 

37°  20 .'5 

121°  38.' 5 

bZ 

1.0 

0.2 

6oooo(v) 

T,F,H 
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(.341) 

(.642) 
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0.8 

2800 

P 
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Kineral 
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40°  20.'7 
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BZ 

1.0 

0.4 

75000 

P 

(.345) 

(.605) 

WHE 

0.8 

2800 

P 

MLC 

Manzanita 

1800 

40°  32.' 2 

121°  33.'7 

LNE 

6.0 

250 

S 

Lake 

(.537) 

(.563) 

ORV 

Qroville 

1180 

39.  33;3 

121°  30.'0 

BZNE 

1.0 

0.75 

100000 

P 

(.555) 

(.500) 

GZNE 

15 

100  3000 

P 

PCC 
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91 

37°  30.'0 

122°  22.'9 

bZ 

1.0 

0.2 

50000(v) 
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(.500) 
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Priest 

1187 

36°  08:5 
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bZ 

1.0 

0.2 
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Valley 

(.142) 

(.665) 

& 
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Paraiso 

363 

36°  19.'9 

121°  22.'2 

bZ 

1.0 

0.2 

7000(V) 

70000(v) 
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Springs 

(.332) 

(.370) 

SAO 

S.  A.  G.  0. 

230 

36°  45  .'9 

121°  26.'7 

bZ 

1.0 

0.2 

100000(V) 

IjFjI'm 

10  km  SW 

(.765) 

(.44^) 

HZ 

J.  15 

-  2x10' (V) 

H,F 

of  Hollister 
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Table  2  (continued) 

« 

MAGNETIC  TAPE  CHANNEL  ASSIGJWENT  -  FROM  JANUARY  14,  1966  TO  JANUARY  23,  1967. 
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5.  TELESEISM  .ANALYSIS 
B.A.  Bolt  and  W.  H,  Bakun 

Many  seismograms  from  teleseisms  are  st’  'ied  in  the  course  of  the 
work  of  the  Seismographic  Station.  The  seismograms  are  used  for  many 
purposes,  some  of  a  routine  nature  (e.g. ,  hypocenter  location  and  tsu¬ 
nami  warning)  and  some  for  spscialized  research.  As  an  illustration  of 
the  latter,  recently  Berkeley  recordings  of  FKP  waves  were  used  in  the 
design  of  an  experiment  involving  the  construction  of  a  sjecial  array 
in  Antarctica  to  study  core  phases.  Two  examples  of  research  at  the 
Seismographic  Station  with  a  general  interest  are  now  presented. 

First,  with  the  important  assistance  of  Dr.  O.W.  Nuttli,  empirical 
P  wave  delay  times,  relative  to  the  Berkeley  station,  have  been  calcu¬ 
lated  for  12  stations  in  central  and  northern  California.  Eight  sta¬ 
tions  of  the  telemetry  network  were  involved.  A  number  of  important 
results  emerged: 

(a)  The  telemetry  stations  must  be  divided  into  two  groups;  the 
northern  group  (CLS,  PRC,  MHC,  SCC,  VIT,  BEK)  evidently  lie  on  a  reason¬ 
ably  uniform  crust  and  upper  mantle  while  the  southern  telemetry  groups 
(LLA,  PRS,  PRl)  axe  situated  on  a  basement  which  significantly  differs 
from  the  northern  one.  The  stations  in  the  Central  Valley  of  California 
show  relative  delays  which  are  quite  unlike  those  of  the  telemetry  sta¬ 
tions.  (Jamestown  has  not  yet  been  studied. ) 

(b)  The  particular  division  of  northern  California  into  regions 
based  on  the  contouring  of  station  residuals  published  by  J.  Cleary  and 
A.L.  Hales  is  probably  not  correct. 

(c)  The  azimuth  anomaly  functions  provide  a  new  way  of  determining 
regional  variations  in  the  Earth ' s  upper  mantle.  The  causal  velocity 
inhomogeneities  must  occur  at  depths  of  order  100  km. 

As  a  separate  study  a  frequency  analysis  of  crustal  reverberations 
of  dilatational  waves  from  telese.isms  is  in  progress;  the  goal  is  to 
obtain  information  on  the  fine  structure  of  the  crust  of  central 
California. 
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A  search  wap  made  of  the  magnetic  tape  library  for  earthquakes 
suitable  to  such  a  study.  From  June  1964  to  October  1966,  over  110 
earthquakes  were  recorded  by  the  Berkeley  array  with  magnitude  ^  5  6111(1 
focal  depth  60  km.  Almost  all  were  in  the  Circum -Pacific  Belt  and  90 
per  cent  had  epicentral  distances  of  55  tc  8o°.  Of  these  earthquakes,  40 
had  magnitude  5-8,  which  was  judged  to  be  the  cutoff  for  earthquakes 
having  a  signal  to  noise  ratio  of  2  to  1.  Eighteen  were  from  the  South 
Pacific  area,  12  were  from  South  America,  and  12  were  from  the  Japan 
Islands-Kurile  Islands  area.  These  earthquakes  were  examined  in  detail 
and  7  were  selected  using  the  following  criteria: 

(1)  Signal  tc  noise  ratio  of  4  to  1. 

(2)  P-phase  duration  cf  10-20  seconds,  \ ncomplicated  by  other 
phase  arrivals  (pep  or  pP\  This  would  yield  10-20  cycles 
of  pure  P-phase  crustal  reverberations. 

Five  of  the  earohquakes  selected  were  from  South  America  and  2  from  Japan. 
(PcP  follows  P  too  closely  at  the  Berkeley  array  for  aU  "  uth  Pacific 
earthquakes  of  interest  in  this  study. )  In  making  uhese  selections,  there 
was  some  indication  that  P-phase  amplitudes  from  earthquakes  having  azi¬ 
muths  relative  to  Berkeley  in  the  range  NT5°W  to  S75°W  (Marianas  Islands- 
Solomon  Islands  region)  were  less  than  expected  from  comparisons  of  South 
American,  Japanese ,  and  other  South  Pac if ic  shocks. 

A  comparison  of  the  traces  recorded  by  stations  across  the  array 
(Figure  l)  indicates  that  the  response  to  a  dilatational  wave  incident 
at  the  base  of  the  crust  below  central  California  is  generally  uniform 
with  moderate  variations  from  station  to  station.  This  suggests  a  gener¬ 
ally  uniform  crustal  structure  under  the  central  California  array.  Of 
special  interest  are  the  phases  al,  a2,  a3,  and  a4  marked  on  the  Mount 
Hamilton  trace  (MHC).  These  phases  can  be  seen  to  some  degree  on  all 
the  traces.  The  deep-focus  character  of  the  source  indicates  that  these 
phases  do  not  originate  in  crustal  reverberations  at  the  source.  The  near 
normal  incidence  (20-25°)  of  the  waves  at  the  base  of  the  crust  suggests 
that  these  prominent  pnases  in  the  vertical  component  are  not  P  to  S  con¬ 
versions  at  the  base  of  the  crust.  One  explanation  of  the  large  amplitude 
of  the  a.L  phase,  especially  on  w>e  GCC,  PRS,  and  FRI  traces,  is  that  the 
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source  consists  of  at  least  tvo  pulses,  separated  by  an  interval  of 
several  seconds.  The  phases  a2  and  a4  could  then  be  interpreted  as 
crustal  reverberations  of  P  and  al. 

Another  possibility  is  that  al,  a2,  a3,  and  a4  are  reflections 
from  a  crustal  interface  of  the  first  P  impulse.  Interference  with  the 
PcP  phase  could  account  for  a3  not  being  prominent.  The  3  to  4  second 
period  would  then  indicate  a  reflecting  surface  5  to  10  kilometers  deep. 
Further  analysis  of  this  and  other  earthquakes  should  limit  the  range  of 
proposed  crustal  models  for  central  California. 
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6.  AFTERSHOCK  SEQUENCES 

T.V.  McEviliy  and  K.B.  Casaday 

A  primary  research  objective  of  the  large-scale  telemetered  array  instal¬ 
lation  and  operation  is  the  study  of  properties  of  aftershock  sequences  in  the 
area.  In  the  original  contract  period,  only  one  significant  sequence  (magni¬ 
tude  5  or  greater)  occurred  and  was  analyzed  by  Udias  (1965). 

During  the  period  covered  by  this  report,  three  aftershock  sequences  have 
been  studied  in  detail.  These  are  the  Corralitos  sequence  of  November,  1964 
(McEviliy,  1966a),  the  Antioch  (Mt.  Diablo)  sequence  of  September,  1965 
(McEviliy  and  Casaday,  1967),  and  the  Parkfield  sequence  of  June-July,  1966 
(McEviliy,  1966b  and  Bakun  et  al. ,  i960). 

Studies  of  aftershock  sequences  are  considered  important  because  they 
reflect  the  nature  of  continued  stress  release  in  the  focal  region  of  larger 
earthquakes  and  contain  valuable  information  on  the  focal  mechanisms  of 
related  earthquakes.  In  general,  the  data  acquired  through  these  studies 
constitute  basic  material  on  sequence  characteristics  which  must  je  taken 
into  account  in  any  theory  of  earthquake  generation. 

When  an  aftershock  sequence  begins  within  the  Berkeley  array,  it  is 
detected  and  located  quickly  by  means  of  the  telemetered  data.  The  procedure 
at  Berkeley  is  then  to  dispatch  mobile  vans  to  the  inferred  epicentral  region 
so  that  hypocenter  determinations  will  be  more  accurate  and  very  small  magni¬ 
tude  events  can  be  monitored.  Coordination  is  also  established  with  the 
various  other  organizations  in  the  area  concerned  with  seismic  studies. 

The  first  sequence  to  be  studied  in  detail  at  the  Berkeley  Seismographic 
Station  was  the  Salinas -Watsonville  sequence  of  1963.  The  methods  applied 
and  the  parameters  established  to  describe  that  sequence  are  similar  to  those 
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measured  In  the  three  more  recent  sequences  and  therefore  a  comparison  of 
four  sequences  can  be  made  below. 

Corralltos  sequence  -  November,  1964.  A  sequence  of  more  than  100  after¬ 
shocks  with  magnitudes  as  low  as  -0.1  was  recorded  following  a  magnitude  5.0 
earthquake  on  November  16,  1964  on  the  San  Andreas  fault  zone  of  central 
California.  The  sequence  was  monitored  in  detail  by  three  temporary  seismo- 
graphic  stations  at  distances  less  than  15  km  and  by  the  surrounding  tele¬ 
metry  array.  Nearly  all  of  the  35  earthquakes  which  could  be  located  clus¬ 
tered  in  a  focal  region  about  4  km  in  diameter  at  a  depth  near  12  km  and 
exhibited  uniform  first  motion  radiation  patterns.  The  slope  of  the  cumula- 
ti”e  frequency  vs.  magnitude  curve  (a  measure  of  large  aftershocks)  was  -0.66. 
First  motion  fault  plane  solutions  are  consistent  with  the  right  lateral 
transcurrent  motion  characteristic  of  the  San  Andreas  fault.  Exceptions  to 
this  uniform  radiation  pattern  in  the  concentrated  focal  region  occurred  near 
the  times  of  two  large  aftershocks  apparently  on  another  fault  about  5  km  away. 

Antioch  sequence  -  September  1965.  A  foreshock-aftershock  sequence 
associated  with  an  earthquake  of  magnitude  4,9  on  September  10,  1965  in  the 
San  Francisco  Bay  region  was  recorded.  The  sequence  was  monitored  by  nearby 
Berkeley  network  stations  and  by  one  temporary  station  in  the  epicentral 
region.  Precise  hypocenters  were  determined  for  29  of  the  shocks  and  these 
clustered  in  a  small  focal  region  with  dimensions  on  the  order  of  a  few 
kilometers  and  at  an  average  depth  of  1?  km.  Wichln  the  focal  region  a 
north-south  trend  was  observed  with  shocks  apparently  migrating  down  and 
then  back  up  the  trend.  First  motion  radiation  patterns  were  uniform 
throughout  the  sequence,  and  a  P-wave  fault  plane  solution  for  the  main 
shock  gave,  as  a  possible  source  a  north-south  trending  fault  with  right- 


lateral  motion.  The  log  cumulative  frequency  o.  shocks  vs.  magnitude  plot 
exhibits  a  slope  of  -O.78. 

Parkfield  sequence  -  June-July,  1966.  Two  earthquakes  of  magnitude  5.1 
and  5-5  Fhook  the  Parkfield  area  in  southern  Monterey  County,  California  on 
June  28,  1966.  Moderate  damage  was  incurred  and  ground  rupture  was  evident. 
These  large  shocks  were  preceded  by  two  foreshocks  of  magnitude  3.1  and  1.3. 
Following  the  second  large  shock,  a  major  aftershock  sequence  began,  with 
half  of  the  shocks  occurring  in  the  first  2k  hours.  On  June  29,  about  39 
hours  later,  a  third  large  shock  of  magnitude  5.0  occurred.  The  sequence  has 
continued  for  several  months  with  some  200  shocks  of  magnitude  1.5  or  greater 
recorded. 

At  the  time  of  the  initial  shocks,  there  were  three  stations  monitoring 
the  area  at  distances  less  than  100  km. ,  Priest  and  Paraiso  of  the  Berkeley 
network  and  Gold  Hill  of  the  U.  S.  Geological  Survey.  As  the  sequence  began, 
the  Berkeley  Seismographic  Station  dispatched  a  mobile  van  to  the  area  and 
assisted  ESSA  personnel  in  installing  some  of  their  mobile  units.  Within  a 
few  days  approximately  20-25  temporary  stations  were  installed  by  various 
groups  within  the  region. 

Although  the  data  acquired  from  this  sequence  are  undergoing  continu¬ 
ing  study,  the  bulk  of  it  has  been  subjected  to  the  standard  analysis  and 
is  reported  here. 

Locations  for  the  foreshocks,  main  »n^cke,  and  for  108  aftershocks,  all 
of  magnitude  *  l.f  have  been  reasonably  well  determined,  and  lie  in  a  region 
some  5  km  in  width  and  35*^0  km  in  length  along  the  San  Andreas  fault  son e 
and  coinciding  with  the  zone  of  surface  cracking  and  maximum  intensity.  The 
average  depth  appears  to  be  8  km  or  less,  although  this  is  subject  to  further 
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revision.  Cumulative  frequency  vs.  magnitude  gives  a  -0.49  slope,  indicat¬ 
ing  a  large  proportion  of  higher  magnitude  aftershocks. 

A  study  of  aftershocks  during  a  short  time  interval  within  the  sequence 
was  carried  out  and  shows  a  similar  spatial  distribution  of  shocks  as  the 
entire  sequence.  Uniform  migration  in  time  or  space  of  the  shocks  was  not 
observed. 

The  surface  cracking  exhibited  right  lateral  displacement  where  measur¬ 
able,  continuing  throughout  the  sequence.  The  extent  of  cracking  is  greater 
than  expected  for  a  magnitude  5*  5  earthquake  (loch  ',  1958)  if  we  assume  that 
the  entire  zone  of  surface  rupture  evident  the  following  morning  (10  hours 
later)  was  involved  in  the  main  shock.  Tocher's  data  would  associate  a  mag¬ 
nitude  7.0  earthquake  with  a  35  to  40  km  fault  rupture.  Intensity  data, 
though  basically  imprecise,  indicate  an  elongated  source  (  Bakun  et.  al. , 
1966).  Acceleration  along  the  surface  rupture  zone  some  30  km  from  the 
main  shock  epicenter  reached  0.5  g  (Cloud,  1966)  along  the  rupture  and  maxi¬ 
mum  values  decrease  away  from  the  rupture  zone.  The  relationohip  of  the 
large  zone  of  surface  breakage  to  the  total  energy  release  in  the  main 
shock  is  not  clear. 

A  P-wave  first  motion  focal  mechanism  analysis  was  applied  to  the  two 
large  earthquakes  of  June  28  using  some  25  stations  at  distances  up  to  several 
hundred  km.  Both  events  had  the  same  radiation  pattern.  The  solution,  with 
no  inconsistent  points  and  observations  in  all  four  quadrants,  gives  a  verti¬ 
cal  strike-slip  right  lateral  fault  strike  of  N33°W,  coincident  with  the 
surface  rupture  and  fault  trace  in  the  area.  Strike  and  dip  are  well  deter¬ 
mined,  to  within  5°  allowable  variation,  by  the  data.  The  motion  appears  to 
have  an  upward  component  on  the  west  side,  at  20°  from  the  horizontal,  with 
the  data  allowing  no  more  than  .’0°  variation  in  this  value. 
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Comparison  01  Aftershock  Studies,  In  the  process  of  studying  central 
California  aftershock  sequences,  attempts  have  been  made  to  develop  parame¬ 
ters  with  which  to  characterize  a  particular  sequence.  The  various  sequences 
studied  can  then  be  compared.  In  the  studies  considered;  the  following  parame¬ 
ters  have  been  employed:  (l)  ^ain  shock  magnitude;  (.2)  relative  magnitude  of 
the  largest  aftershock  as  compared  to  the  main  shock;  (3)  slope  of  the  cumu¬ 
lative  frequency  vs.  magnitude  curve  (also  a  measure  of  proportion  of  large 
aftershocks);  (4)  average  focal  depth  of  the  aftershocks;  (5)  long  dimension 
of  the  focal  region  (defined  as  the  region  of  aftershock  activity).  Tabulated 
results  for  t’..'-  Salinas-Uatsonville,  Corralitos,  Antioch,  and  Parkfield 
sequences  are  shown  below. 


Relative  Magni¬ 

Long  Dimension 

Main  Shock 

tude  of  largest 

Ave. 

of  Focal 

Sequence 

Magnitude 

Aftershock 

Slope 

Depth 

Region 

Watsonville 

5.4 

.9 

-.41 

7 

30 

1963 

Corralitos 

5.0 

.7 

-.  66 

12 

4 

1964 

Antioch 

4.9 

.6 

-.78 

12 

6 

1964 

4 

Parkfield 

5.5 

.9 

-.49 

8 

35-40 

1 

1966 

S§ 

From  these 

measures  of 

the  aftershock  sequences  it 

appears 

that  in  all 

characteristics 

used,  Corralitos  and  Antioch 

are  similar  and  Salinas- 

3 

Watsonville  and 

Parkfield  are  similrr.  The  geographic 

location  of  the  sequence 

n 

s 


appears  immaterial.  Salinas -Watsonville  and  Corralitos  are  in  nearly  the  same 
area.  Antioch  iB  in  a  much  less  seismically  active  area  than  the  remaining 
three. 

It  seems  from  these  parameters  that  a  larger  focal  region  and  greater 
number  of  large  aftershocks  occur  together,  not  surprisingly  in  conjunction 
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with  the  larger  min  shocks.  The  magnitude  difference  of  one-helf  unit  between 
similar  pairs,  however,  seems  somewhat  small  relative  to  the  order  of  associated 
variations  in  the  other  sequence  characteristics.  Sequences  associated  with 
larger  main  shock  magnitudes  are  also  characterized  by  shallower  average  depths. 

Whether  or  not  these  results  reveal  two  distinct  types  of  earthquake 
sequences  in  central  California,  two  points  on  a  continuous  variation  in 
parameters  with  magnitude  >r  chance  observations  from  a  wide  variety  of 
sequence  behaviour,  will  have  to  be  determined  on  the  basis  of  subsequent 
aftershock  studies.  Complementing  such  further  studies  with  observations 
of  ultramicroearthquake  activity,  strain  variations,  S-vave  characteristics, 
and  spectral  information  will  add  valuable  data  in  the  search  for  understand¬ 
ing  the  mechanisms  of  an  earthquake. 
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7.  CHASE  V  AND  THE  CHICO  EARTHQUAKE 
C.  Lomnitz 


The  original  purpose  of  the  Chase  V  experiment  was  as  follows: 

(1)  Calibration  of  travel -time  curves  used  to  locate  earthquakes  in 
the  offshore  Cape  Mendocino  area.  The  seismicity  of  this  region  is  of 
interest  in  relation  to  the  structure  at  the  intersection  of  the  San 
Andreas  fault  with  the  Mendocino  Escarpment. 

(2)  Determination  of  crustal  structure  and  seismic  velocities  in 
northern  California. 

In  1965,  a  seismologies!,  experiment  involving  the  detonation  of  at 
least  500  tons  of  explosives  in  the  Pacific  Ocean  was  proposed  by  B.A. 

3olt  and  C.  Lomnitz  to  the  Office  of  Naval  Research.  The  coordinates 
suggested  and  subsequently  approved  by  ONR  were  40°  36'  N,  125°  45'  W,  in 
the  center  of  an  active  seismic  region  about  75  miles  west  of  Cape  Mendocino. 

Operation  Chase  V.  An  explosion  of  surplus  ammunition  of  approxi¬ 
mately  1000  tone  cf  TNT  equivalent  occurred  on  May  24,  1966  off  the  coast 
of  California.  The  explosives  had  been  loaded  into  an  expendable  Liberty 
ship  and  towed  to  a  location  off  Cape  Mendocino.  Because  of  operational 
problems,  the  detonation  occurred  earlier  and  at  a  different  location  than 
originally  planned. 

We  undertook  the  task  of  planning  and  coordinating  the  seismic  experi¬ 
ment  on  shore,  with  the  participation  of  recording  units  from  the  California 
Department  of  Water  Resources,  ESSA,  Stanford  Research  Institute,  U.  S. 
Geological.  Survey,  University  of  British  Columbia,  University  of  Nevada,  and 
University  of  Oregon.  Observations  obtained  by  ARPA,  the  California  Institute 
of  Technology,  and  Columbia  University  have  also  been  utilized. 

The  sites  of  the  temporary  stations  were  selected  to  provide  long  pro¬ 
files  oriented  on  the  main  axis  of  *be  U.  C.  Berkeley  seismic  array,  plus 
two  shorter  profiles  terminating  in  the  Oroville  and  Jamestown  stations. 

Since  the  detonation  occurred  about  130  km  to  the  south  of  the  intended 
shotpoint,  the  geometry  of  the  experiment  was  substantially  altered 
(Figure  l). 
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Two  hours  before  the  CHASE  V  explosion  a  natural  earthquake  of  magni¬ 
tude  4.5  occurred  in  the  vicinity  of  Chico,  California,  This  earthquake  was  well 
enough  located  (Figure  l)  to  provide  comparisons  of  travel -times  and 
amplitudes  with  Chase  V;  accurate  epicentral  location  was  facilitated  by 
the  temporary  stat  3ns  which  had  recorded  both  events. 

Table  1  gives  a  list  of  selected  readings  for  Chase  V  and  for  the 
Cuico  earthquake.  The  focal  coordinates  and  Berkeley  magnitudes  are  also 
given.  This  list  wf  3  combined  with  reported  arrival  times  from  other  sta¬ 
tions  to  produce  the  travel- time  graph  for  P  shown  in  Figure  2  over  the 
range  of  distances  from  0°  to  10°. 

The  following  preliminary  conclusions  are  suggested  by  these  results: 

(1)  There  are  systematic  differences  between  the  seismic  /elocities 
computed  from  the  separate  sets  of  data.  In  general,  these  ..ifferences 
can  be  attributed  to  the  circumstance  that  seismic  paths  from  Chase  V  were 
largely  off  shore  or  grazing  the  coast  line,  while  the  paths  for  the  Chico 
earthquake  were  all  continental. 

(2)  The  average  apparent  surface  velocity  of  Pn  for  Chase  V  was  7. 97 
km/sec  (Figure  2)  which  is  close  to  the  value  proposed  many  years  ago  by 
Byerly  for  central  California, 

(3)  The  average  apparent  surface  velocity  of  P  for  the  Chico  earth¬ 
quake  was  significantly  lo^er,  7.58  km/sec  (Figure  2). 

(4)  Seismic  paths  which  cross  the  roots  of  the  Sierra  Nevada  are  sig¬ 
nificantly  delayed  in  both  events;  thus  the  times  to  Tinemaha,  Mina  (Nevada), 
China  Lake,  Isabella,  and  Eureka  are  late  by  an  average  2.6  seconds  with 
reference  to  either  average  Pr  travel  time  curve,  of  7.97  km/sec  for  Chase 

V,  and  7.58  km/sec  for  the  Chico  earthquake  (Figure  2). 

(5)  Paths  which  are  purely  or  mostly  offshore  tend  to  produce  earlier 
arrivals,  which  appear  to  be  consistent  with  a  P^  velocity  of  the  order  of 
8.2  km/sec. 

It  is  also  of  interest  to  compare  the  character  of  the  signals  for  the 
two  events.  Figure  3  shows  a  playback  from  the  magnetic  tape  system  using 
two  array  stations  which  were  located  at  comparable  distances  from  the  source 
(about  400  km).  In  each  case,  the  upper  trace  (a)  is  the  unfiltered  signal; 
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the  other  traces  were  obtained  by  narrow  band -pass  filtering  at  frequencies 
of  (b)  4  cps,  (c)  1.8  cps,  and  (d)  1  cps. 

For  Chase  V,  a  large  proportion  of  the  energy  is  contained  in  the  1.8 
cps  band,  which  corresponds  to  the  frequency  of  the  bubble  pulse.  No  sig¬ 
nificant  S  phase  is  apparent.  It  is  also  of  interest  that  later  phases  in 
the  P-group,  such  as  P*  and  P,  are  much  more  pror  Lnent  in  the  Chico  earth¬ 
quake  than  in  Chase  V. 

The  large -amplitude  wave  train  observed  in  the  late  portion  of  the 
Chase  V  record  presumably  corresponds  to  a  surface  wave  which  propagates 
in  the  ocean  bottom  sediments,  since  its  velocity  (about  1.75  km/sec)  is 
higher  than  expected  lor  a  T  phase  (  1.45  km/sec}. 

The  structure  of  the  initial  portion  of  the  seismograms  for  both 
events  may  be  compared  in  Figure  4,  a  high-speed  playback  with  the  same 
settings  as  for  Figure  3. 

A  more  complete  discussion  of  these  results  will  be  given  at  the 
forthcoming  meeting  of  the  Seismological  Society  of  America  in  Santa 
Barbara  in  a  paper  now  under  preparation  by  C.  Lomnitz  and  B.A.  Bolt. 
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Table  1 


Chase  V 

(Suirpliei  location) 


Chico  Earthquake 

( Preliminary  location) 


Latitude 

39°  28'  N 

39°  43'  N 

Longitude 

125°  48'  W 

121°  44'  w 

Depth 

3,750  ft 

25  km 

Water  depth 

12,500  ft 

- 

Bubble  pulse  period 

0.576  sec 

- 

Magnitude 

4.7 

4.5 

Origin  time 

O5.49-O6.8 

03-49-54. 3 

P  Arrival  Times 

0B3  Ocean  Bottom 

O5.49-31.2 

- 

PAR  Pt.  Arena 

34.2 

03-50-23.4 

ARC  Areata 

37.4 

29.5 

UKI  Ukiah 

38.1 

17.0 

FLA  Plantation 

39.7 

22.1 

P3K  Paskenta 

75.1 

06.8 

^LS  Calistoga 

45.6 

14.8 

PRC  Point  Reyes 

47.8 

24.2 

RUM  Rumsey 

48.8 

11.9 

SHS  Shasta 

51. 5 

11.3 

QUE  Quebec 

- 

19.3 

TAN  Tango 

- 

21.0 

KIL  Kilo 

- 

21.5 

BRK  Berkeley 

54.8 

24.2 

BKS  Berkeley  (Byerly) 

55.0 

24.4 

CM?  Concord 

5<5.4 

22.5 

ROM  Romeo 

- 

22.5 

JUL  Juliet 

- 

25.2 

ORV  Oroville 

57.0 

02.0 

PCC  Pilorcitos  Creek 

56.6 

30.0 

MIN  Mineral 

57.7 

05.8 

GCC  Granite  Creek 

50-04.0 

36.1 

MHC  tit.  Knailton 

04.8 

32.2 

(continued) 
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Table  1  -  continued 

P  Arrival  Times 

Chase  V  Chico  Earthquake 


KFO  Klamath  Falls 

10.  T 

35.5 

KRC  Karris  Ranch 

11.7 

40.9 

SPR  Spreckels 

11.4 

42.8 

JAS  Jamestown 

12.9 

28.5 

LLA  Llanada 

16.8 

43.9 

PRI  Priest 

23.9 

51.2 

FRE  Fresno 

28.7 

47.5 

COR  Corvallis 

28.8 

51-11.0 

EUR  Eureka,  Nevada 

58.5 

02.0 

BMO  Blue  Mountains  Obs. 

51-08.6 

- 

TIN  Tinemaha 

50-43.5 

50-59.7 

WDY  Woody 

43.4 

51-04.7 

ISA  Isabell , 

49.0 

09.7 

SBC  Santa  Barbara 

45.2 

17.8 

FTC  Ft.  Tejon 

51.4 

16 

CLC  China  Lake 

54.8 

15.2 

PAS  Pasadena 

51-04 

28 

GLC  Goldstone 

05.7 

24 

MWC  Mt.  Wilson 

04.6 

29 

RIV  Riverside 

- 

35 

FLC  Palomar 

20.7 

47 

BRT  Barrett 

^0.0 

56 

VIC  Victoria,  B.C. 

23.3 

- 
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8.  SEI3CCITY  OF  THE  REGION 
B.A.  Bolt,  C.  Lomnitz  and  T.V.  McEvilly 

We  have  taken  the  point  of  view  that  one  of  the  major  research 
objectives  for  which  the  Central  California  Seismic  Array  could  be  used 

is  the  study  of  seismicity  patterns  in  the  Coast  Ranges  of  California. 

« 

The  recent  locations  c2  array  stations  have  been  selected  in  large  part 
to  give  optimum  control  in  hypocenter  determinations  in  the  seismically 
active  regions. 

Figure  1  shows  a  composite  map  of  earthquake  epicenters  of  magni¬ 
tude  2.5  and  greater  in  the  central  California  region  (aftershocks  not 
included),  as  determined  by  stations  of  the  Berkeley  network  during  the 
years  1962,  1963*  and  1964  (first  three  quarters).  Results  for  1965  are 
now  being  processed.  The  following  points  of  special  interest  deserve 
to  be  mentioned.. 

(1)  In  the  Coast  Ranges  most  of  the  seismic  activity  during  1962, 
1963,  1964  has  occurred  in  two  main  regions  both  within  the  Gabilan 
Range:  (a)  the  Salinas -San  Juan  Bautista -Watsonville  region  west  of 
Hollister;  (b)  the  Pinnacles -Bear  Valley  region  south  of  Hollister. 

Both  regions  are  situated  to  the  west  of  the  San  Andreas  fault 
trace.  We  believe  that  this  is  probably  a  genuine  offset  but  the  effect 
of  systematic  errors  in  epicenter  determinations,  due  to  different  veloci¬ 
ties  east  and  west  of  the  fault,  is  beii  investigated. 

(2)  North  of  Corral it os,  no  seismicity  is  apparent  (magnitude 
*2.5)  on  the  San  Andreas  fault  until  the  fault  goes  out  to  sea  off 

the  San  Francisco  Peninsula.  A  possible  exception  is  an  earthquake  of 
magnitude  4.6  in  the  Butano  Ridge  area,  about  15  km  west  of  the  fault 
trace.  The  recorded  seismicity  off  the  San  Francisco  Peninsula  is  inter¬ 
preted  as  a  remanent  effect  of  the  aftershock  sequence  following  the  1957 
Paly  City  earthquake  (magnitude  5-7). 

(3)  No  earthquakes  of  magnitude  greater  than  2.5  were  recorded  on 
the  San  Andreas  fault  between  Marin  County  end  Pt.  Arena  with  the  possible 
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exception  of  the  single  event  southeast  of  Pt,  Arena,  some  13  km  east 
of  the  fault  trace, 

(4)  In  the  Diablo  Range  near  the  Junction  of  the  Hayward  and 
Calaveras  faults  there  are  active  seismic  loci  tfhich,  in  this  case, 

are  very  closely  correlated  with  the  surface  strike  of  the  major  faults. 
The  seismic  activity  continues  northwest  into  the  area  of  Calaveras 
Reservoir,  north  of  San  Jose. 

(5)  Earthquakes  along  the  Hayward  fault  tend  to  cluster  slightly  to 
the  west  of  the  fault  trace.  Most  of  the  activity  has  been  in  Alameda 
County. 

(6)  Activity  along  the  Calaveras  fault  was  concentrated  principally 
between  Mt,  Diablo  and  Livermore,  slightly  east  of  the  fault  trace.  Thie 
corresponds  to  the  epicentral  area  of  the  Antioch  earthquake  of  September 

10,  1965. 

(7)  There  is  an  indication  of  sporadic  seismicity  on  the  eastern 

flank  of  the  Coast  Ranges.  These  epice’-'  .*•  mostly  around  Panoche 

Valley  and  into  the  western  San  Joaquin  Valley;  this  seismicity  deserves 
close  attention  because  of  the  potential  hazard  to  the  California 
Aqueduct. 

(8)  There  are  scattered  epicenters  in  the  northern  Coast  Ranges, 
north  of  Carquinez  Straits,  particularly  in  Sonoma,  Lake,  and  Rapa 
Counties.  Geologic  mapping  of  faults  is  difficult  in  this  area  due  to 
lack  of  surface  contacts. 

(9)  Determination  of  focal  depths  of  earthquakes  plotted  in  the 
Coast  Ranges  indicates  in  this  interval  of  study  (1962-1964)  at  least 
80  per  cent  of  the  better -determined  shocks  are  no  more  than  about  5  km 
deep.  Focal  depths  in  excess  of  10  km  are  exceptional. 
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9.  PROPERTIES  OF  S -WAVES  AT  NEAR  DISTANCES 
J.  Filson 


An  effort  has  been  made  to  study  the  elastic  shear  wave  from  sources  at 
near  distances.  The  research  has  been  directed  toward  determining  the  ground 
particle  motion  due  to  the  arrival  of  the  shear  wave  at  these  distances;  and 
toward  determining  the  mechanism  of  propagation  which  will  produce  the 
observed  motion.  The  latter,  of  course,  involves  properties  of  the  medium 
and  the  source. 

To  60lve  the  observational  problem,  the  responses  of  two  horizontal  and 
one  vertical  Press-Ewing  seismometers,  located  at  Berkeley,  are  recorded  on 
magnetic  tape.  These  seismometers  have  a  free  period  of  30  seconds  and  like 
damping;  they  have  been  calibrated  and  their  response  characteristics  are 
considered  identical.  Upon  playback,  the  horizontal  components  are  referred 
to  a  resolving  device  which  electronically  rotates  the  seismometer  response 
to  the  longitudinal  and  transverse  directions  with  respect  to  the  station 
to  epicenter  line.  Thereupon,  the  vertical,  longitudinal,  and  transverse 
signal  components  are  passed  through  operational  amplifiers  whose  vesponsf 
characteristics  are  such  that  the  response  of  the  seisn  ometer-amplifier  sys¬ 
tem  is  rendered  constant  with  respect  to  displacement  between  30  and  .1 
seconds.  Particle  motion  diagrams  in  the  vertical-transverse,  vertical- 
longitudinal,  and  transverse -longitudinal  (horizontal)  planes  are  then  con¬ 
structed  from  recorder  traces.  In  some  cases  additional  filtering  is  required 
to  suppress  the  microseisms.  Examples  of  ground  particle  motion  due  to  the 
arrival  of  the  shear  wave  from  sources  within  40  km  are  shown  in  Figure  (l). 

In  consideration  of  a  theory  which  might  predict  the  observed  particle 


motions,  an  isotropic  homogeneous  half  space  propagating  plane  harmonic  waves 
serves  as  a  simple  model.  For  further  simplification,  the  theoretical  S 
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motion  is  resolved  into  vertical  SV  and  horizontal  SH  components.  At  the 
point  of  observation  on  the  surface,  the  SH  motion  is  reflected  without  change 
in  amplitude  while  at  certain  angles  of  incidence  a  portion  of  the  SV  motion 
may  be  converted  into  pure  longitudinal  motion,  a  P  wave  being  generated. 

Upon  this  conversion  at  the  surface,  the  angle  between  the  P  wave  front  normal 
and  the  surface  normal  is  greater  than  for  the  reflected  SV  wave.  ThuB  there 
is  a  critical  angle  of  incident  SV  above  which  there  is  no  generated  P  wave. 

The  surface  particle  motion  for  angles  of  incidence  less  than  this  critical 
angle  s  linear;  however,  for  incidence  angles  greater  than  this  critical 
angle,  the  SV  component,  changes  phase  with  respect  to  the  SH  component  and 
the  resulting  particle  motion  describes  a  three  dimensional  figure  in  Bpace. 

A  program  was  written  for  the  plane  wave,  half -space  model,  to  compute 
and  draw  the  theoretical  particle  motion  for  any  angle  of  incidence  and  for 
any  ratio  of  original  SH  and  SV  components.  Figure  (2),  drawn  by  the  program, 
shows  the  three  planes  of  motion  as  the  angle  of  incidence  passes  through  the 
critical  angle  (here  about  37°)  in  5°  steps.  Notice  that  the  motion  is  linear 
up  to  35°,  then  elliptical  at  40°,  linear  again  at  45°  with  no  longitudinal 
component,  and  finally  elliptical  again  at  50°,  For  all  angles  of  incidence 
greater  than  50*  the  motion  is  non-linear.  The  P  and  S  wave  velocities  used 
in  drawing  Figure  (2)  are  8.1  and  4.6  km/sec,  respectively. 

Now  note  the  particle  motion  of  the  S  wave  of  earthquake  #4  in  Figure  1. 

It  is  nearly  linear  with  a  definite  longitudinal  component.  For  this  linearity 
to  be  explained  by  the  model,  the  focal  depth  would  have  to  be  greater  than  ths 
epicentral  distance  of  37  km.  This  earthquake,  which  is  one  of  the  Antioch 
sequence  discussed  in  Section  6,  is  probably  not  deeper  than  12  km.  Thus  this 
direct  S  wave  should  give  rise  to  non-linear  motion  at  the  surface.  In  this 


case,  a  shallow  surface  layer  of  greatly  reduced  shear  velocity,  say  1.2  km/sec, 
is  required  to  yield  the  observed  linear  motion.  If  this  zone  exists  its  thick¬ 
ness  would  have  to  be  such  that  it  would  have  little  effect  on  the  longer  period 
S  phases  from  earthquakes  at  intermediate  distances.  These  waves  give  rise  to 
elliptical  SV  motion  as  reported  by  Nuttli  and  Whitmore  (1962).  Although  the 
remaining  two  examples  in  Figure  1  appear  complicated,  they  and  most  of  the  S 
arrivals  from  epicenters  within  80  km  are  more  linear  than  not,  and  certainly 
not  elliptical  as  predicted  by  the  plane  wave,  half  space  model.  Thus  this 
model  is  considered  to  fail,  except  for  the  possible  variation  employing  the 
surficial  low  velocity  zone. 

The  model  presently  being  studied  is  one  of  a  buried  impulsive  source  in 
a  single  layer,  over  a  half  space.  Again,  the  half  space  and  the  layer  are 
isotropic  and  homogeneous;  the  layer  modeling  the  crust,  the  half  space,  the 
mantle.  Here  one  must  consider  the  direct,  reflected,  and  conical  waves  and 
the  P  wave  when  generated  at  the  upper  and  lower  boundaries  of  the  layer. 
Although  theories  of  handling  this  problem  have  been  developed  in  general  by 
Brelchovskikh  (i960),  Cagniard  (1962),  and  Heelan  (1953)*  the  specific  case 
the  surface  motion  due  to  the  arrival  of  the  conical  S  wave  remains  to  ''e 
worked.  It  is  toward  the  solution  of  this  problem  that  this  research  is  now 
directed. 
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10.  CRUSTAL  STRUCTURE  ESTIMATION  WITHIN  THE  ARRAY 

T.V.  McEvilly 

The  University  of  California  seismic  array,  located  between  the  west¬ 
ern  continental  margin  and  the  Sierra  Nevada  range,  clearly  spans  a  zone 
where  complex  transitions  in  crust  and  upper  mantle  structure  should  occur. 

In  addition,  the  array  is  divided  b  the  San  Andreas  fault  zone  through  the 
Coast  Ranges,  with  the  resultant  lateral  discontinuity  of  granitic  rocks  to 
the  west  and  the  Franciscan  series  to  the  east.  It  is  not  surprising  that 
anomalous  features  of  wave  propagation  are  discovered  in  studies  utilizing 
wave  arrivals  at  stations  in  this  area.  Understanding  and  considering  these 
anomalies,  however,  is  of  prime  importance  in  fully  utilizing  the  capabili¬ 
ties  of  the  large-scale  array  in  seismic  research;  for  example,  in  obtaining 
maximum  accuracy  of  epicenter  location  both  for  local  and  distant  earthquakes. 
On  the  other  hand,  the  advantages  of  array  recording  can  be  brought  to  bear 
on  the  problem  cf  resolving  the  underlying  structural  variations  producing 
the  anomalous  propagation. 

Seismic  studies  of  several  types  have  produced  strong  c-vidence  of  non- 
uniform  crust  and/or  upper  mantle  structure  in  the  array  area.  Evernden 
(1953,  1951*,  1955)  observed  variations  from  great  circle  paths  in  directions 
of  approach  of  Rayleigh  waves  recorded  in  central  California.  Otsuka  (1966a, 
b),  working  with  the  array,  found  cyclic  anomalies  in  apparent  velocity  and 
direction  of  approach  for  teleseismic  P-waves  which,  :.n  general,  were  indica¬ 
tive  of  either  very  rapid  crustal  thickening  eastward  or  a  low  velocity  zone 
in  the  upper  mantle  thickening  eastward.  Bolt  and  Nuttli  (1966)  1  simi¬ 
lar  variations,  relative  to  Berkeley  observations,  in  teleseismic  P  wave  resid¬ 
ua’  at  individual  stations  of  the  array  in  northern  California.  Refraction 
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studies  in  central  California  also  suggest  that  crust-mantle  configuration 
is  complex.  Interpretations  have  suffered  from  a  lack  of  good  control  on  the 
branch  of  th'  travel-time  curve.  Studies  by  Byerly  (1939)>  Cameron  (1961), 
Ryall  (1962),  Mealy  (1963),  Eaton  (1963),  and  Hamilton  et  al.  (1964),  have  pro- 


km/sec  difference  across  the  San  Andreas  zone,  velocities  from  7.6  ‘o  8.0 
km/sec,  and  crustal  thicknesses  of  20  to  30  km. 

During  the  contract  period  summarized  in  this  report  an  effort  wa6  made 
to  improve  the  apparent  uncertainty  in  crustal-u  per  mantle  structural  inter¬ 
pretations  in  central  California  Coast  Ranges.  The  method  devised  utilizes  the 
seismic  array,  implemented  by  some  fifty  temporary  stations  that  have  been 
placed  in  operation  over  the  past  fo’a  yearn  in  various  studies,  as  a  tri¬ 
dimensional  refraction  array  for  sources  consisting  of  about  twenty  accurately 
timed  blasts  at  three  locations  and  precisely  located  hypocenters  in  four 
aftershock  sequences  (see  smother  section  in  this  report).  The  data  are  a 
suite  of  travel  times  between  more  than  100  station-source  pairs.  Distances 
range  up  to  nearly  300  km. 

The  Interpretation  method  is  a  least  squares  procedure  of  minimizing 
residuals  in  the  observed  travel  tim_-s  by  systematically  varying  structural 
parameters.  Crustal  velocity  is  approximated  by  a  linear  increase  with  depth. 
Surface  velocity,  gradient,  crusted  thickness,  upper  mantle  velocity,  Mobo  dip 
and  curvature  are  variable  parameters.  Preliminary  results  using  only  blast 
data  are  summarized  by  McEvilly  (1966).  This  in*  Mai  work  emphasized  the  dif¬ 
ficulty  of  recording  P  arrivals  from  explosions  in  the  area.  P  arrivals, 

n  g 

however,  gave  strong  evidence  for  a  change  of  about  0.6  km/sec  in  shallow 
crustal  velocity  across  the  San  Andreas  zone,  with  higher  velocities  on  the 
southwestern  side.  Combining  all  data  for  upper  mantle  structure  yielded  a 
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depth  of  about  22  km  to  the  Moho  in  the  center  of  the  array  with  appreciable 
east  dip  and  a  Pr  velocity  of  7.83  km/sec.  This  result  suffers  from  the  poor 
quality  of  Pn  arrivals  from  explosions.  The  study  has  been  expanded  to  include 
large  aftershocks  in  sequences  studied  in  detail  using  temporary  stations  at 
very  short  distances  (~ 10  km)  so  that  errors  in  focal  depths  and  origin  times 
are  not  significant.  The  required  program  modifications  for  variable  source 
depth  have  been  made  and  the  additions]  data  are  now  being  readied  for 
analysis. 

Overall  general  conclusions  can  be  drawn  from  the  various  studies 
related  to  crust  and  upper  mantle  structure  in  central  California.  There  is 
strong  evidence,  as  well  as  geophysical  reason,  for  general  inland  dip  on  the 
Moho.  Unfortunately,  quantitative  determination  of  dip  magnitude  is  difficult 
since  the  effect  of  the  lateral  velocity  decrease  across  tie  San  Andreas  zone, 
when  viewed  through  an  averaging  process  using  the  entire  a^ray,  is  virtually 
equivalent  to  that  produced  by  crustal  thickening  inland,  or  even  seaward  dip 
on  the  top  of  a  low  velocity  layer  within  the  mantle.  The  simplest  manner  of 
resolving  the  ambiguities  calls  for  determining  the  shallow  structural  proper¬ 
ties  and  subsequently  "stripping"  these  known  effects  to  investigate  deeper 
features.  The  indications  of  a  lateral  velocity  decrease  of  the  order  of  0.5 
km/6ec  in  the  shallow  crust  northeast  of  the  San  Andreas  are  fairly  good  and 
the  details  can  be  mapped  by  refraction  studies,  P^  arrivals,  when  kno^r. 
crustal  velocities  are  used,  will  contribute  information  on  depth  and  config¬ 
uration  of  the  Moho.  Present  interpretations  favor  a  20-  to  25-km  crust  with 
several  degrees  of  dip  inland  and  possible  crustal  thickening  to  the  southeast. 
Better  Pr  arrivals  within  the  array  are  required  to  clarify  this  structure. 


The  velocity  distribution  in  the  upper  mantle  must  be  inferred  from  the  tele¬ 
seism  azimuth  and  velocity  anomalies,  taking  into  account  the  structure  above 
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At  the  present  stage  of  interpretation,  it  is  clear  only  that  the  seaward 
dipping  velocity  inversion  in  the  upper  mantle,  when  considered  in  conjunction 
with  a  moderately  inland  dipping  Moho,  produces  velocity  and  azimuth  anomalies 
similar  to  these  observed.  In  terms  of  anomaly  the  dipping  low  velocity  lid 
is  essentially  equivalent  to  a  lateral  inland  decrease  of  about  0.5  km/sec  in 
velocity  in  the  upper  mantle. 

In  summary,  within  the  area  of  the  continental  margin  spanned  by 
the  array,  it  appears  that  velocities  decrease  landward  in  both  the  crust 
(where  t;.e  jump  appears  associated  with  the  San  Andreas  fault  zone)  and  the 
mantle.  Total  change  across  the  array  may  be  as  much  as  0.5  km/sec  in  each 
case.  The  Moho  appears  to  exhibit  a  few  degrees  of  landward  dip  and  lie  at 
an  average  depth  of  20  to  25  km  within  the  array.  Systematic  studies  in  com¬ 
bining  refraction  data,  teleseism  anomalies,  and  possibly  surface  wave  dis¬ 
persion  are  required  for  further  definition. 
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